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ABSTRACT 

We examine the effects of the engulfment of planets by giant stars on the 
evolution of late-type stars. We show that the rate at which dynamo-generated 
magnetic energy is being released exceeds 10% of the wind kinetic energy when 
the orbital angular momentum of the engulfed planet is more than ten times 
the angular momentum of the star as it leaves the main sequence. A significant 
enhancement in the mass-loss rate may be expected in this case, due to the 
formation of cool magnetic spots. We use the existing sample of extrasolar planets 
to estimate that at least 3.5% of the evolved solar-type stars will be significantly 
affected by the presence of planetary companions. 



Subject headings: planetary systems - stars: evolution - stars: mass-loss - stars: 
magnetic fields - stars: statistics - stars: low mass, brown dwarfs 
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1. Introduction 

Radial velocity surveys of nearby stars have been enormously successful in discovering 
extrasolar planets. These surveys have demonstrated that a significant fraction of the stars 
in the solar neighborhood have massive planets with orbital radii that are substantially 
smaller than Jupiter's (e.g. Marcy & Butler 2000; Udry et al. 2001; Butler 2001; Vogt et al. 
2002). The engulfing of these planets by their parent stars is an inevitable outcome of 
stellar evolution (either directly or due to tidal interaction) in the red giant branch (RGB) 
or asymptotic giant branch (AGB) phase (e.g. Livio & Soker 1984; Siess & Livio 1999a,b). 
The capture of planetary-mass companions has been suggested to be the potential cause 
for the high rotational velocities (V^otsini > 10 km s^^) observed in some field red giants 
(Stefanik et al. 2001) and blue horizontal branch stars (e.g. Peterson et al. 1983; Soker & 
Harpaz 2000). 

In the present work we examine specifically the effects of planet/brown dwarf "swal- 
lowing" on the rotation rates of and mass loss from giants. In particular, we estimate the 
expected statistics of influenced stars. 



2. Rotation and Mass Loss in Evolved Late- Type Stars 

When a star in the RGB or AGB phase engulfs a planet/brown dwarf the latter deposits 
its angular momentum into the giant's envelope. The resulting rotation frequency, as a 
fraction of the critical stellar rotation frequency oOk = {GM^/ R^)^ (where M*, R^, are the 
stellar mass and radius respectively) is given by 



u / Mr \ f M„,^\ ^ f kl\ (a 



Here Mc is the companion mass, M^nw is the mass of the giant's envelope, MkgR"^ is the star's 
moment of inertia and a is the initial separation between the giant and the companion. Since 
stars like the Sun rotate only at a fraction of a percent of their critical rate (on the main 
sequence) and this fraction is further reduced (by conservation of angular momentum) to 

(-) <0.l(-) f— , (2) 

as the stars ascend the RGB or AGB, we see from eq. (1) that even Jupiter-mass planets can 
have a very significant effect on the rotation rate of late-type stars. Even when the increase 
in central condensation is taken into account, typically, {cu/cUkjAGB ~ (0.01 — 0.1){u!/u!k)MS 
(e.g. Eriguchi et al. 1992). 
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The effects on mass loss rates are more difficult to estimate because they involve the 
processes of generation of magnetic activity (especially in spots) and mass loss via radiation 
pressure on molecules and dust, both of which are rather poorly understood. Nevertheless, 
it has been shown semi-quantitatively that the spiraling-in process of the companion inside 
the giant's envelope can result in a dynamo generation of a magnetic field (Regos & Tout 
1995). 

Here we shall adopt a simplified, more heuristic approach (similar to that of Soker & 
Harpaz 1992) to estimate at what rotation rates we can expect significant effects on mass 
loss. Dynamo-gcncratcd magnetic fields arc limited by the requirement that the amplification 
timescale be of the order of the buoyant rise time of magnetic flux tubes (e.g. Parker 1979, 
§8.7). In a typical stellar aui dynamo the amplification timescale is given approximately by 
(e.g. Zeldovich, Ruzmaikin & Sokoloff 1983) 



(3) 



where a = iuJo/3 is determined by convection [£ is the mixing length and ujo the surface 
angular velocity), k is the wave number and Va; is a measure of differential rotation. The 
critical wave number kc below which dynamo waves are amplified is given by {i^r is the 
turbulent viscosity) 

fcf ~ ai?,Vcj/2z/2 . (4) 

Adopting the same values as Zeldovich ct al. (1983), Wuj = uJo/{0.3R^), = iv^/S, where 
Vc is the convective velocity, and taking the wave number to be equal to the critical one, we 
find 



1 



For typical AGB star parameters, Vc ^ 10 km s~^, <£>~ -R*/4, the obtained average 
amplification time is of the order of Tq — 9rrot, where Trot is the rotation period. The local 
rise time of a fiux tube is given approximately by (Parker 1979) 

Avc AirpAvc 

where Va is the Alfven speed, B is the magnetic field strength, A r/2 (where r is the 
radial distance from the center; e.g. Soker & Harpaz 1999) is the local scaleheight and p is 
the density. The density in the convective zone of an AGB star can be well approximated 
by (Soker 1992) 

Menv \ f -R* V V r _3 
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Equating the flux tube rise time to the regeneration and amphfication time Ta yields a 
value for the magnetic field 

B . 1.8 (^] ^ (^] ^ (^) ^ ( ,] ^ (^] ^ (^] G . 

\MqJ \300RqJ \0.1MqJ VlOkms-V V IQ-'y VlOO R©/ 

(8) 

Magnetic energy is thus being released at a rate 
1 f 

Eb ^ — — Anr'^dr 

Ta J StT 

^ 1028 (M^\ ( Vc \ ^ (^^\ (MA ' ' erg s^ (9) 

The kinetic energy associated with the mass loss from an AGB star is of order 

where M is the mass-loss rate and V^j is the wind velocity. An examination of eqs. (9) 
and (10) reveals that for the companion to have an appreciable effect on mass loss (i.e. to 
have Eb > 0.1£^wind) it must increase the angular momentum of the envelope (and thereby 
u!o/u!k) by an order of magnitude (or more). Consequently, in our estimate of the fraction of 
stars that are affected significantly by planetary companions, we shall adopt as the defining 
condition 

Jp > lOJ* , (11) 

where Jp is the orbital angular momentum of the planet, and J* is the angular momentum 
of the star as it leaves the main sequence [eq. (2)]. Note that this is probably a reasonable 
criterion for the planet having a significant effect on stellar evolution even considering the 
uncertain nature of the effects on mass loss. In fact, given that the magnetic field in stellar 
spots can be ~ 10^ times stronger than the average field (e.g. Priest 1987), condition (11) 
may be too stringent, since even a field lower by an order of magnitude than that given by 
eq. (8) will result in spots that are significantly cooler than their surroundings (magnetic 
pressure larger than the photospheric pressure). 



3. Statistics on Importance of Planets 



We define the following quantities 

u = 
H = 



log{a/AU) 
\og{M/Mj) 



(12) 
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where a is the orbital separation, M is the minimum planetary mass, and Mj is Jupiter's 
mass. In Fig. 1 we show the known 75 planets (taken in January 2002 from the Extrasolar 
Planets Encyclopaedia compiled by J. Schneider, at http: / /www. obspm.fr/encycl/cncycl. html) 
in the u-fi plane. In planetary systems where more than one planet is known wc mark the 
most massive planet by an asterisk, while the others are marked by the + sign. Several 
authors have attempted recently to obtain unbiased distributions of planets as a function 
of u and jj, (e.g. Lineweaver & Grether 2002; Mazeh & Zucker 2002; Tabachnik & Tremaine 
2002; Armitage et al. 2002). In view of the uncertainties, we adopt the following simple form 
for the fraction of systems dN within a box di/d/i (for the parameter space indicated below) 

dN = C {I + au){l- (5 ij)dndu , (13) 

for —1.5 < < 0.5 and ji < 0.5i/ + 0.75. Here C, a and (3 are parameters. This form 
is the same as the one used by Lineweaver & Grether (2002), who found (their figures 3 
and 4) q; ~ 0.8 and (5 ~ 0.8. Using the well populated parallelogram (i.e. excluding areas in 
the V — IX plane with large gaps; see also Zucker & Mazeh 2002) marked by lines 1 to 4 on 
Figure 1, we estimate a ~ 0.3 and j3 ^ 0.5. Examining in addition the results of Mazeh & 
Zucker (2002), Tabachnik & Tremaine (2002), and Armitage et al. (2002), and fitting where 
necessary their functional forms to ours, gives average values of <a>~ 0.6 and </5>~ 0.5, 
which we shall now adopt. 

The value of C can be determined for fixed values of a and (3 in the following way. From 
a sample of 1,200 stars discussed by Vogt et al. (2002; where more details can be found) 
44 planets and 4 brown dwarfs have been detected. Hence 3.7% of the surveyed stars were 
found to have planets. The 75 known extrasolar planets (as of the end of January 2002) 
are distributed in 67 planetary systems. Using this ratio we find that the fraction of stars 
being detected to harbor planetary systems is ~ 3.3%. If we examine the well populated 
area of Figure 1 (bounded by the four straight fines v — 0.5; ji — 0.5i/ + 0.75; v — —1.5; and 
IX = O.Biy), we find that out of the 67 planetary systems, 48 are inside this parallelogram. 
This corresponds to ~ 2.4% of the stars having their massive planet within this region. 
Considering that the right side of this region is close to the detection limit (no planets are 
found near u ^ 0.4 and fi ~ 0.3), wc take the detection fraction within this region to be 
Ng — 2.5% of all stars. For stars with a metallicity above solar, [Fe/H] > 0, the detection 
fraction is ~ 2.6 times higher than that for the entire sample (Vogt et al. 2002). Hence, in 
the same scaling-region the fraction is Ng ~ 6.5%. The large increase with metallicity in 
the probabifity of harboring a planet was studied by Gonzalez (1997) and Reid (2002). The 
latter author argues that most stars with a metallicity [Fe/H] > 0.3 have planets around 
them. We therefore take Ng — 0.025 for all stars and Ng — 0.065 for metal rich stars, and 
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find the value of C from the expression 



/0.5 r 
du 
-1.5 -'0. 



0.5I/+0.75 



c 

32 



Ns^C 



dii{l + au){l- Pii) 



(48 - 24q; - 6/3 - 5aP) = 0.9094a (14) 



where the last equality was obtained by substituting the average values a = 0.6 and (3 — 0.5. 

In the previous section we derived as our condition for a significant effect on mass-loss 
from the parent star, Jp > 10 J*. We noted, however, that it is possible that even a lower 
angular momentum might play some role, especially since the star loses angular momentum 
as it evolves (Sokcr 2001). Soker further argues that planets of masses as low as O.OlMj are 
sufficient to influence the mass loss geometry from AGB stars. An indirect hint that even 
a relatively slow rotation may influence the mass loss rate comes from the distribution of 
stars on the horizontal branch in globular clusters (GCs). D'Cruz et al. (1996) show that 
the distribution of stars on the horizontal branch in the HR diagram can be explained by 
variations in the mass loss rate on the RGB, by a factor which depends on the metallicity 
(ranging from 1.5 in metal rich GGs to > 5 in metal poor GGs). We should note that red 
horizontal branch stars also rotate (Stefanik et al. 2001). However, using the results of Behr 
et al. (2000) for M13, Sokcr & Harpaz (2000) showed that the average angular momentum 
of the progenitors of hotter horizontal branch stars is larger than that of cooler ones. It 
is thought that rotation is the main factor that influences the mass loss rate variations 
(R. Rood, private communication, 2002), with the exception of very high mass loss rates 
which may result from interactions with stellar binary companions. Since main sequence 
stars in GCs are expected to rotate very slowly (since they have experienced a long period of 
angular momentum losses due to magnetic activity), if the rotation conjecture is correct, then 
very slow rotation is sufficient to influence mass loss on the RGB. Our condition, therefore, 
may be somewhat conservative. 

The present angular momentum of the Sun is 1.7x lO^^g cm^ s~^. A planet with a Jupiter 
mass at 1 AU from a solar type star has ~ 50 times the present angular momentum of the Sun. 
Taking as atypical value J* = Jq, our condition in equation (11) reads {M/Mj){a/ AU)^/^ > 
0.2. For a more general form Jp > r}J^ — tjJq, our condition reads 



In addition, wc limit the planet to a < 2 AU (or i/ < 0.3), in order for tidal forces to be 
able to bring it into the RGB or AGB star's envelope. These conditions are represented (for 
7] = 10) as the two thick lines in Fig. 1. Planets to the left of and above these lines affect the 
evolution of their parent stars signiflcantly (neglecting eccentric orbits of planets at larger 
orbital separations). 



/x + 0.5z/>log(VlO)-0.7 = j, 



(15) 
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Not including the few systems above the hne marked '2' in Fig. 1, the total fraction of 
planetary systems that are likely to significantly influence the evolution of the central star 
at late stages is therefore, 



Note that this integral is meanigful for j„i < —0.7. For jm > —0.7 the thick diagonal line in 
Fig. 1 crosses the upper limit as given by line '2.' Performing this integral for a = 0.6 and 
P — 0.5, and substituting for C from equation (14), we find 



Equation (17) shows that for our strong condition of = 10, we find that more than 
3.5% of all solar-type stars are significantly affected (in terms of rotation rate and mass 
loss) by planets or brown dwarf companions. This number increases to more than 9% for 
metal rich (above solar) stars. Relaxing the condition to = 1, i.e. when the planet's orbital 
angular momentum equals the present one for the sun, increases the numbers to be > 8% and 
> 21% respectively. All of these numbers should probably be regarded as lower hmits, since 
we have neglected the effects of planets at 2 AU< a <5 AU in eccentric orbits (which could 
also be pulled into the giant's envelope by tidal interaction; see also Debes & Sigurdsson 
2002), and used only the minimum masses (as well as not having included the planets above 
line '2' in Fig. 1). 



We have shown that the swallowing of planetary or brown dwarf companions can have 
significant effects in the late stages of the evolution of late-type stars. Specifically, we 
considered rotation and mass loss. On the basis of the properties of the observed extrasolar 
planets, we estimated that ~ 4-10% of the stars should experience enhanced mass loss due 
to the effects of planets. 

We acknowledge useful comments by an anonymous referee. N. S. was supported by a 
Celerity Foundation Distinguished Visiting Scholar grant at the University of Virginia, and 
by a grant from the US-Israel Binational Science Foundation. He also asknowledges the 
hospitality of ST Scl. 
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N{total) = Ns{OA - 1.2jm + O.SfJ . 



(17) 



4. 



SUMMARY 
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Fig. 1. — 75 known extrasolar planets are presented in the minimum-mass-semimajor-axis 
logarithmic plane. In planetary systems where more than one massive planet has been found, 
the most massive one is marked by *, while the others by +. The four hues marked 1 to 4 
define the area used for our scaling of fraction of planets. The two thick lines represent our 
condition for a planet to significantly influence the mass loss rate of its parent star (see text) . 



